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ABSTRACT 

Fast Radio Bursts (FRBs), are millisecond radio signals that exhibit dispersion larger 
than what the Galactic electron density can account for. We have conducted a 1446 
hour survey for Fast Radio Bursts (FRBs) at 145 MHz, covering a total of 4193 
sq. deg on the sky. We used the UK station of the LOFAR radio telescope - the 
Rawlings Array - , accompanied for a majority of the time by the LOFAR station 
at Nangay, observing the same fields at the same frequency. Our real-time search 
backend, ARTEMIS, utilizes graphics processing units to search for pulses with dis¬ 
persion measures up to 320 cm -3 pc. Previous derived FRB rates from surveys around 
1.4 GHz, and favoured FRB interpretations, motivated this survey, despite all previ¬ 
ous detections occurring at higher dispersion measures. We detected no new FRBs 
above a signal-to-noise threshold of 10, leading to the most stringent upper limit yet 
on the FRB event rate at these frequencies: 29 sky -1 day -1 for 5 ms-duration pulses 
above 62 Jy. The non-detection could be due to scatter-broadening, limitations on the 
volume and time searched, or the shape of FRB flux density spectra. Assuming the 
latter and that FRBs are standard candles, the non-detection is compatible with the 
published FRB sky rate, if their spectra follow a power law with frequency (oc u a ), 
with a > +0.1, demonstrating a marked difference from pulsar spectra. Our results 
suggest that surveys at higher frequencies, including the low frequency component of 
the Square Kilometre Array, will have better chances to detect, estimate rates and 
understand the origin and properties of FRBs. 

Key words: instrumentation: miscellaneous — pulsars: general 
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1 INTRODUCTION 

Higli-time resolution radio astronomy is entering a new era 
due to two developments. Firstly, new radio telescopes are 
being built that combine very high instantaneous sensitivity 
with wide fields of view, as required for blind searches of 
short duration events of transient radio emission. Secondly, 
high performance computing (HPC) enables us to process 
the data extremely quickly, providing the opportunity to 
discover and react to such events, enabling localization and 
classification. 

Our quest to explore the dynamic nature of the radio 
sky is underpinned by the realization that bright, short dura¬ 
tion radio bursts must originate from extreme physical pro¬ 
cesses taking place around compact sources. Neutron stars 
offer prime examples of such phenomena, being reliable emit¬ 
ters of short duration (typically 50 ms or less) radio pulses 
( Lorimer fe Kramer||2005 1. These signals are observed once 
per rotational period for the majority of ordinary pulsars, 
but can be as infrequent as a few pulses per day in the case 


of Rotating Radio Transients (RRATS; McLaughlin et al. 


2006). In addition to neutron stars, a new class of extraordi¬ 


nary sources has been identified in high time resolution radio 
astronomy which highlights the importance of real-time pro¬ 
cessing. The first of these sources - that have now come to 
be known as Fast Radio Bursts (FRBs) - was discovered 


by Lorimer et al. (2007), who reported the detection of a 


burst using the Parkes radio-telescope, with evidence sup¬ 
porting large, extragalactic distance. Following this discov- 


2015). We stress that all these events except Petroff et 


Thornton et al. 

|2013| Spitler et al. 2014; Burke-Spolaor &; 

Bannister 2014; 

Petroff et al.|2015; Ravi, Shannon, & Jame- 


al. (2015) were discovered via offline processing of recorded 


data, therefore no immediate follow-up observations to lo¬ 
calize them were possible. The origin of FRBs is unclear. 


Proposed explanations include flaring magnetars (Popov & 
Postnov 2013), binary neutron star mergers (Tot ani| 2013[), 
gravitational collapse of neutron stars to black holes (Falcke 
& Rezzolla 20141, emissions from companions of extragalac- 


tic pulsars ( |Mottez fe Zarka 2014), and at least one Galac¬ 
tic proposition in the form of nearby flare stars (Loeb et al. 


20141. These associations can be strengthened or weakened 


by comparing the rates of these events to the observed rates 
of FRBs. Current best limits from the observed FRBs are 
given by 


Thornton et al. 
volumetric rate of 10 -3 


(20131 as 10 4 

T" 


sky day , or a 


gal * year 

In general, an FRB will be dispersed by the free electron 
content of the medium it propagates through. As a result, 
the lower frequency components of a broadband burst arrive 
later than their high frequency counterparts, the delay at a 
given frequency v being proportional to the line integral of 
the electron column density along the line of sight (known 
as the dispersion measure) and v~ 2 . There is some evidence 
that the dispersion measure (DM) of one of the aforemen¬ 


tioned FRBs - the Keane et al. (20121 burst - may be ex 


plained by the Galactic distribution of electrons (Bannister 
|fe Madsen|2014[ ). If FRBs are indeed extragalactic in origin, 
they can be used to quantify the ionized matter in the inter- 
galactic medium (IGM), thereby allowing us to determine 
the baryon content of the universe. This could help solve 
the long-standing ‘missing baryon problem’ in cosmology, 


where there is a discrepancy between the observed and the 
expected quantities of baryons (Persic k, Salucci||1992 Me- 


Quinn|[2014 1. Answering this question is one of the science 


objectives of the Square Kilometre Array (SKA Macquart 
et al.|2015 ). 


There is an expectation that astrophysical bursts will 
emit across a broad range of radio frequencies. The observed 
signal-to-noise ratio of an FRB can be maximized by correct¬ 
ing for the dispersion induced delays across the frequency 
band of a telescope, in a well established process from pul¬ 
sar astronomy called dedispersion. In the specific context 
of a blind search for FRBs, where the DM is not known, 
the main challenge is to perform the incoherent dedisper¬ 
sion transform on the data. By this, we mean the process of 
producing from time series of total power in a large num¬ 
ber of narrow frequency channels, time series of total power 
over a large range of DM values, revealing potential FRB 
candidates. This is a computationally expensive task, and 
performing it in real time requires specialized software run¬ 
ning on high performance computers with accelerated hard¬ 
ware. In the last few years, implementations of the accel¬ 
erated dedispersion transform that are capable of real-time 


processing of typical data streams have emerged (Armour 
et al.|20T2 Barsdell et al.|2012| Magro et al.||20lT l and are 

gradually being applied to surveys. 

Comprehensive arguments in favour of searching for in¬ 
dividual bright and dispersed pulses have been present in 
the literature and applied in practice for many years (see, 
for example, McLaughlin & Cordes 2003). Yet despite the 


discoveries mentioned above, we are still in the early stages 
of putting together a complete picture of the dynamic na¬ 
ture of the radio sky. The most up-to-date limits hitherto on 
fast transient sources at 150 MHz can be found in ICoenenl 
|et al.| ( p014[ ), who use the results of the LOFAR Pilot Pulsar 
Survey (LPPS) to derive an upper limit on the FRB rate of 
150 sky -1 day -1 . All aforementioned FRBs have been dis¬ 
covered at a frequency around 1400 MHz. Upper limits to 
the FRB rate at different frequencies provide constraints on 
the luminosity distribution and the spectral characteristics. 
Assuming the flux density spectra of FRBs follow a power 
law, then 


S(v) oc , 


(1) 


where S(u) is the flux density of the burst as a function of 
frequency u, a is the spectral index. Low frequency FRB sur¬ 
veys are also being carried out at the Murchison Widefield 
Array (MWA), a low frequency dipole array located in West¬ 
ern Australia operating between 80 and 300 MHz. |Trott,| 
Tingay, & Wayth (2013) suggest several expected FRBs per 


observing week, providing a clear framework for estimating 
FRB rates depending on their intrinsic spectral index and 
scattering properties. There is currently no published detec¬ 
tion of an FRB from the MWA in the literature. 

The Low Frequency Array, LOFAR, is currently the 
most sensitive radio telescope in the world at frequencies 
between 30 and 250 MHz. It is a prototypical example of a 
‘next-generation radio telescope’ and a pathfinder towards 
the Square Kilometre Array (SKA). The telescope is de¬ 
scribed in detail in van Haarlem et al. (20131. Discovery of 


pulsars and fast transients is one of the key science goals of 
LOFAR, as outlined in |Stap pers et al. (2011). We present 
here details and first results of a programme to use one inter- 
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national LOFAR. station to search for FRBs, in an effort to 
characterize the dynamic nature of the radio sky at LOFAR 
frequencies. In this paper, we describe the system developed 
to conduct fast transient surveys and the results of the first 
survey. 

This paper is organized as follows. In ^2j we describe 
the project, including details of the hardware and software, 
and the signal processing algorithms, and in $3 we describe 
our survey for fast transients using ARTEMIS, before con¬ 
cluding in f|4] 


2 THE ARTEMIS PROJECT 

We have put together the hardware and software that is re¬ 
quired to continuously monitor the large LOFAR. fields of 
view and detect FRBs in real time using general purpose 
computing on graphics processing units (GPUs). The scien¬ 
tific potential and the technical challenge of real-time FRB 
discovery were the prime motivating factors behind the Ad¬ 
vanced Radio Transient Event Monitor and Identification 
System (ARTEMIS) project, with an emphasis on applying 
HPC techniques to data from next generation telescopes. 
The primary goals of the ARTEMIS project therefore con¬ 
sist of: 

• Enabling the science goals of the SKA and its pathfind¬ 
ers in the area of surveys for new pulsars and exotic fast 
transients. 

• Exploiting many-core technologies to minimize capital 
costs of hardware, whilst delivering energy-efficient solutions 
for streaming processing in radio astronomy. 


2.1 Technical description 


We use the high-band antennas (HBA) of a single interna¬ 
tional LOFAR station, which consists of 96 dual-polarization 
antennas covering a frequency range between 110 and 
250 MHz. The ARTEMIS installation currently uses four 
12-core servers, each equipped with a single Fermi NVIDIA 
GPU card located on site at the Rawlings Array in the LiK 
(a similar setup exists at the international LOFAR station 
in Nan§ay, France and an ARTEMIS server for offline pul¬ 
sar data processing exists in Jiilich, Germany). The servers 
receive data through a broadband (10 Gbps) switch, which 
is also responsible for sending the data back to the LOFAR 
correlator in the Netherlands, for normal operations of the 
International LOFAR Telescope (ILT). Each LOFAR station 
generates a total of 3.2 Gbps of beamformed data, which 
corresponds to a sky bandwidth of approximately 48 MHz 
using 16-bit sampling. The station beamformer produces a 
total of 244 beamlets each with the frequency and time res¬ 
olution specified in Table [T] Estimated values of the system 
equivalent noise are also given, for an observation with a 
total bandwidth of 6 MHz and time resolution of 327.68 ps, 
by scaling the sensitivity values, as given in Table B.3. of 


van Haarlem et al. (20131 for a source at 30° declination, to 


the aforementioned bandwidth and sample time. The sensi¬ 
tivity value quoted for 150 MHz is in good agreement with 
the value from the LPPS survey ( Coenen et al.|2014 l using 
uncalibrated LOFAR stations, scaled to this survey. 

For this survey, groups of 30 and 31 beamlets are 


Table 1. Specifications of an international LOFAR station 


Frequency 

Beamlet FWHM 

Sensitivity 

(MHz) 

(deg) (Jy, Ai/=6 MHz, 5t= 327.68 ps) 

30 

9.9 

1566 

60 

4.0 

1079 

120 

2.5 

33 

150 

2 

27 


Elevation range 

30° to 90° 

Maximum number of beamlets 

244 

Maximum total sky bandwidth 

47.65625 MHz 


Sampling rate 

200 or 160 Msamples/s 

Beamlet frequency resolution 

195.3125 or 156.25 kHz 

Beamlet time resolution 

5.12 or 6.4 ps 


serial pipeline 



Figure 1. A block diagram of the ARTEMIS software, outlining 
the main processing modules, buffers and asynchronous tasks. 


pointed in a single direction on the sky, covering contigu¬ 
ous frequency channels to form a total of 8 broadband 
(~ 6 MHz) beams. These can be formed within a circle of 
~10 times the beamlet FWHM, due to the analog tile beam¬ 
forming of LOFAR. HBA arrays described in Stappers et al. I 

poTTj ). 

Each ARTEMIS server receives and processes 2 out of 
the 8 formed beams, corresponding to a quarter of the total 
bandwidth. In addition to FRB detection, real-time process¬ 
ing can be used to reduce the data rate from 400 MB/s to 
manageable rates both for storage and further processing, 
useful for other science goals. 


2.2 Software overview and framework 

ARTEMIS provides an extendible and parallelisable soft¬ 
ware backend, which can process telescope data in real time, 
mainly for detecting FRBs. The software, called AMPP, 
is sufficiently generic to serve as an all-purpose telescope 
backend, which is built upon a C++ framework (PELI- 
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capQ, providing a fully configurable client-server architec¬ 
ture. The software pipelines consist of a number of C++ 
modules and configuration files for processing live telescope 
data. The modular design allows data processing tasks to be 
implemented independently of other considerations such as 
data acquisition, transport protocols, data output, and sys¬ 
tem deployment. The modules currently available include 
an implementation of a polyphase filter for raw (coherent) 
telescope data, a Stokes-generating and integration module 
which converts the raw data to Stokes parameters, a radio¬ 
frequency interference (RFI, unwanted and often man-made 
signals) rejection module which operates on incoherent data, 
an incoherent dedispersion transform module, and stream¬ 
ing modules that record data to disk and send the processed 
data via the network for further processing or recording. Fig¬ 
ure [l] shows a schematic of AMPP. As shown, dedispersion 
is an asynchronous task which occurs only when sufficient 
data have been accumulated in a buffer. This task is also 
the most expensive in terms of computation, and therefore 
is accelerated using GPUs. AMPP uses the AstroAcceler- 
ate dedispersion transform algorithm (Armour et al.j2012|, 
which has been favourably benchmarked against other ex¬ 


tremely fast implementations (Barsdell et al. 2012 Magro 
et al.|2011 1. 


2.3 Spectral processing and RFI excision 


In aperture arrays such as LOFAR, beamforming can require 
that spectral channels are formed in the raw digital data 
| Mol fc Romein|2011| ). This is achieved prior to beamform¬ 
ing via a polyphase filter followed by a fast Fourier trans¬ 
form, which produces a number of raw-data frequency sub¬ 
bands. The channelisation requirements of the beamformer 
are generally different to those of the dedispersion trans¬ 
form, therefore an additional channelisation step is necessary 
within the software. Incoherent dedispersion cannot remove 
intra-channel dispersion. As a result, the channel bandwidth 
should not exceed a certain limit, dictated by the maximum 
DM of a search and the minimum duration of the FRBs be¬ 
ing searched for. Determining this limit is standard practice 
in pulsar and fast transient searches (e.g., Keith et al. 2010). 

We have developed a multi-threaded polyphase filter- 
bank module for AMPP, which efficiently generates 2 n spec¬ 
tral channels for each input subband. The number of output 
channels, the windowing function, the number of taps for 
the filter and the number of processing threads dedicated 
to the module are configurable options. In our searches so 
far, we use a filter that provides frequency channels with 
a bandwidth of 3.0518 kHz, by channelizing each incoming 
raw data subband into 64 channels. This increases the sam¬ 
ple time from the original 5.12 ps to 327.68 ps. 

Within the AMPP pipeline, the raw data are passed 
to a module which generates Stokes parameters from the 
incoming complex polarisations, according to: 


I = XX* + YY* 
U = 2Re(Xy*), 


Q = XX* - YY* 

V = —21m (XY*), (2) 


where * refers to the complex conjugate. In the disper¬ 
sion searches for FRBs described here, only Stokes-7 (total 


1 http://www.oerc.ox.ac.uk/~ska/pelican/ 


power) is used. The software however provides the possibil¬ 
ity to buffer all four Stokes parameters. 

RFI is known to yield large numbers of false positives 
in pulsar searches | |Eatough~e t al. 2009). The dedispersion 
transform integrates total power over frequency, thereby re¬ 
moving all frequency resolved information. This imposes the 
requirement to clean the data of RFI before dedispersion. 
Examples of RFI excision techniques for FRB searches have 
been described in Hogden et al. (2012). For our survey, re¬ 
moving RFI in real-time imposes specific constraints both 
on the algorithm (only information from the past data is 
available) and on the processing. The adaptive threshold¬ 
ing algorithm that we developed for AMPP performs the 
following steps: 


(i) Prior to any observation, the spectral bandpass is 
recorded using narrow frequency channels. Individual bright 
channels are rejected through an iterative rejection process, 
and the clean bandpass is modelled using a low-order poly¬ 
nomial. The root mean square (RMS) of the bandpass noise 
is recorded. 

(ii) The search commences. 

(iii) Within AMPP, the total power data are processed 
by the RFI module one time sample at a time, consisting 
typically of a large number of spectral channels. 

(iv) The median of each time sample spectrum is com¬ 
puted, which is insensitive to outliers, and the bandpass 
model is adjusted to that level (the zeroth order coefficient 
is set to the value of the spectrum mean). 

(v) The power in each frequency channel is compared to 
the model at that frequency. If it exceeds a threshold T c h, it 
is given a small weight (in current implementations we set 
the weight to zero). T c h is typically set to 6 times the RMS 
of the current bandpass noise. 

(vi) The mean and RMS of this ‘cleaned’ spectrum are 
stored in a circular history buffer, which can hold values 
extended up to a duration much greater than the longest 
FRB searched for. 

(vii) A running average of the history buffer is computed 
both for the model mean and the RMS, ensuring that the 
model remains up-to-date when it is compared to the next 
observed spectrum, without being affected by short term 
variability. 

(viii) If the median value of the incoming spectrum re¬ 
veals a large disagreement between the model and the spec¬ 
trum, i.e., a jump has occurred in the incoming spectra due 
to a broadband signal, the whole spectrum is assigned a low 
weight. The threshold T sp is typically also set to 6 times the 
calculated RMS of the median. 

(ix) If the number of consecutive spectra that do not 
agree with the model is equal to the size of the history buffer, 
the algorithm accepts that the bandpass model is no longer 
relevant to the data, and sets the mean and RMS to the 
values of the last observed spectrum. 

(x) Finally, for the FRB search, the bandpass model is 
removed from each flagged spectrum. AMPP has a config¬ 
uration setting that allows each spectrum to be scaled to 
maintain a uniform RMS in the post-filtered data, which is 
very useful for the dedispersion transform. 

Figure [2] shows the effect of the RFI excision algorithm 
on 1.3 s of filterbank data from the Rawlings Array. The al¬ 
gorithm described above effectively flags broadband signals 
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Figure 2. Filterbank data from the Rawlings Array. The central images show total power versus time and frequency. The bottom panel 
shows a time series of the total power integrated across frequency, and the panel to the right shows the bandpass averaged over time. In 
the bandpass, the series of parallel horizontal lines mark the boundaries of the LOFAR subbands, which have been further channelised in 
AMPP. The figure on the left shows data before the RFI excision module, including narrowband interference at ~145 MHz which flares 
and sweeps across a narrow frequency range around time bin 680, and a broadband signal around time bin 200. The figure on the right 
shows the same data, post RFI filter. The broadband signal at time bin 200 has been flagged and in this case replaced by noise. The 
narrowband interference has been completely flagged out, revealing two dispersed bright pulses from PSR B0329—54. originally buried 
in the noise. 


which raise the mean of the spectrum. In this example, the 
flagged data are replaced by Gaussian noise matching the 
bandpass model used for the comparison. However, for the 
real-time FR.B searches, we set the weight of these data to 
zero, and replace them in the input to the dedispersion mod¬ 
ule with noise from a chi-squared distribution with four de¬ 
grees of freedom. A narrowband (and in this case frequency- 
swept) signal of interference is also flagged, revealing two 
dispersed pulses of PSR B0329+54 in the processed filter- 
bank data on the right. The RFI excision module can also 
be set to remove the mean of every spectrum, effectively 
removing signals at 0 DM (Eatough et al. 20091. 


channel dispersion. Cordes & McLaughlin (20031 provide a 


2.4 Dedispersion transform considerations 

The dedispersion module used in AMPP was developed as 


part of the Astro-Accelerate library (Armour et al. 2012). 


This module performs the dedispersion transform between 
total power versus frequency and time, and total power ver¬ 
sus dispersion measure and time. The software framework 
provides the infrastructure for setting up tasks on avail¬ 
able GPUs, and GPU dedispersion kernels are used during 
AMPP processing. The details of the dedispersion process 
are beyond the scope of this paper and form the subject 
of a separate publication. Together with the RFI-cleaned 
data, the mean and standard deviation of the baseline noise 
are also propagated through the software pipeline to facil¬ 
itate the detection of statistically significant signals in the 
de-dispersed data. 

Incoherent dedispersion is performed in steps of DM. 
For technical reasons, pertaining to the dedispersion algo¬ 
rithm on a GPU, we maintain a fixed <5DM step across the 
searched DM range. The dedispersion transform in our sur¬ 
vey was therefore carried out at 5DM steps of 0.1 cm -3 pc. 
This discretization results in some loss of signal due to intra¬ 


framework for estimating the signal loss given the DM step 
size, the centre frequency and bandwidth of the observations 
and the width of a Gaussian pulse. Using their equations 12 
and 13 


where 


WDM) = ^ rVrf(i 


c = 6.91 x 10“ 3 c5DM, AumHz 


Wobsi'o, 


(3) 


(4) 


and <5DM = 0.1 cm -3 pc, we estimate a loss in sensitivity 
of less than 15% for 2 ms wide pulses, which drops to 2.5% 
for a 5 ms pulse. 

Regardless of prior beliefs on the origin of FRBs, at a 
sky frequency of 145 MHz, interstellar scattering is a central 


consideration in the design of a search. Bhat et al. (2004) 


compute an empirical relationship between the measured 
scattering timescale r and the DM of a large number of 
pulsars. There is at least an order of magnitude uncertainty 
in the DM-r relationship. Scattering results in broadening 
of the intrinsic FRB width, resulting in a smaller loss of sen¬ 
sitivity, as per eqs. [3] and [4] To investigate relevant values 
of W 0 bs we assume the worst case scenario as regards dedis¬ 
persion (i.e., the least possible scattering per given DM). 
We therefore set r m i n to be an order of magnitude smaller 
than what is given by Bhat et al. (2004). In principle, the 
observed width of a radio pulse is given by the convolution 
of r m i n with the intrinsic width of the burst, although for 
simplicity, we derive the observed width as 


Wobs = Wfnt 


I 

+ 7"n] 


(5) 


where the intrinsic width Wi n t is defined in some reasonable 
way (e.g., FWHM). Another factor that comes into play 
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Figure 3. The efficiency f or a i> — 145 MHz and 

Au = 6 MHz search for FRBs with intrinsic width of 2 and 5 ms, 
given as a function of DM. The curve is computed u sing Eqs. [3| 
|4| |5| a nd setting r m ; n to be 10% of what is given by Bh at et al.| 
(2004). Regardless of the intrinsic FRB width, Wob s is likely dom¬ 
inated by the scattering timescale r m ; n above a certain DM, thus 
rendering the efficiency high even for low intrinsic widths. 


specifically for signals from cosmological distances is cosmo¬ 
logical time dilation, which stretches the pulse by a factor 
of (I + 2 ). However, as discussed in the following sections, 
2 is small in the case of the ARTEMIS survey, so this effect 
can be safely neglected. 

I 11 the LOFAR HBA band at 145 MHz, there is a reason¬ 
able expectation that Wobs will be dominated by scattering 
above a certain DM. Figure [ 3 ] shows the efficiency 
for a DM resolution <5DM=0.1 cm -3 pc, for FRBs with an 
intrinsic width of 2 and 5 ms, based on the arguments above. 
The curves are characterised by a flat segment at the lowest 
DMs, which indicates the DM range for which Galactic scat¬ 
tering is unlikely to substantially broaden the pulse width 
(Wobs £=* Wint). The figure suggests that a DM resolution of 
5DM = 0.1 cm -3 pc maintains high efficiency even at low 
DM values. 

Considering that FRBs typically last a few milliseconds 
sets an upper limit to the DM value searched using incoher¬ 
ent dedispersion. The reason is that intra-channel dispersion 
becomes large with increasing DM value at these low fre¬ 
quencies. Reducing the channel bandwidth is not an option, 
as it worsens the temporal resolution. For the channel band¬ 
width of 3.0518 kHz, the DM smearing at 145 MHz is 8.3/rs 
per unit of DM. We have therefore chosen the maximum 
DM of our search to be 320 cm -3 pc, corresponding to a 
maximum of 2.6 ms of intra-channel smearing. To search for 
higher DM FRBs requires a step of coherent dedispersion, 
which would leave the temporal resolution unaffected. 

The ARTEMIS hardware is capable of processing at 
least 3200 DM values per data stream in real-time, which is 
sufficient given the DM resolution arguments above. In fact, 
this sampling of the DM space was also chosen to demon¬ 
strate the viability of the hardware and software, operating 
at near maximum load. Dedispersed timeseries at neighbour¬ 
ing DM values will show some degree of correlation, however 
the threshold used for detection (10 times the noise RMS, see 


the following section), precludes spurious detections caused 
by this effect. 

Once the filterbank data have been converted to a time 
series per DM, these are decimated to 2, 4, 8, 16, 32 and 
64 times the original sampling time, and, along with the 
original, searched in real time for significant peaks by ap¬ 
plying a 5-sample median filter to preclude the detection 
of single, high S/N peaks. We are currently investigating 
adaptive filtering algorithms for computational and detec¬ 
tion performance. 


3 FRB SEARCHES USING THE RAWLINGS 
ARRAY 


We have currently concluded 1522 hours of drift-scan obser¬ 
vations using the Rawlings Array. These observations were 
accompaniedfor the majority of the time by the French LO¬ 
FAR international station at Nangay. The purpose of the 
second station was to perform a coincidence test on poten¬ 
tial detected FRBs and rule out the possibility that they are 
caused by RFI. 

The array is configured to point 8 beams with a 
bandwidth of 6.055 MHz (for odd-numbered beams) and 
5.859 MHz (for even-numbered beams), to fixed positions in 
the sky on the local meridian, at declinations of 7.9°, 11°, 
13.3°, 15.85°, 19°, 22°, 24.9°and 28.66°, corresponding to 
declinations of known bright pulsars which we use as sanity 
checks. Each beam is about 1.9°wide, and the total amount 
of instantaneous sky coverage is 24 sq. deg. Our sensitivity 
is given by 


S min = 24.7 Jy ( -^= 


( 6 ) 


where 24.7 Jy is the noise RMS in each sample averaged 
across the band. This is computed by correcting the value 
given in Table [l] by the average elevation of our beams 
(58°). Furthermore, 10 is the threshold S/N ratio that we 
used, and D is the decimation factor, compared to the na¬ 
tive 327.68 ps time resolution. This makes us sensitive to a 
peak flux density > 62 Jy for 5 ms-duration pulses. This cor¬ 
responds to a fluence, i.e. energy contained within a pulse, 
of 310 Jy ms. Note that, given the other decimation factors 
we have searched for, our survey has multiple fluence limits, 
ranging from 81 Jy ms for D=1 (pulse durations of 327ps) 
and 650 Jy ms for D=64 (pulse durations of 21 ms). 

The survey searched for FRBs up to a DM of 320 cm -3 
pc, covering a volume of 1.15 x 10 8 Mpc 3 , which was com¬ 
puted in the following manner: Using the NE2001 model of 


electron density in the Galaxy (Cordes & Lazio 2002) and 


making the simple assumption that the potential host galaxy 
contributes a constant DM of 100 cnU 3 pc, we can estimate 
the excess DM due to the IGM in each of our pointings. 
Figure [4] shows the pointings of the survey in Galactic co¬ 
ordinates. Pointings that did not probe outside the Galaxy 
(in the Galactic plane, indicated by red markers in Figure[4| 
were excluded from the volume calculation. Then, following 
|Lorimer et al.| (|2007| ) , we used the results of |Ioka| ( |2003| ) 
and Inoue (20041 to estimate the redshift corresponding to 


the IGM component of the DM. There are considerable un¬ 
certainties in the relationship between redshift and DM. Ac¬ 
knowledging this, we proceed by adopting DM=1200 2 cm -3 
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Figure 4. The pointings of the eight beams of the drift-scan sur¬ 
vey. The red markers represent locations where the Galactic plane 
obscures extragalactic sources, while the green markers represent 
pointings that probed outside the Galaxy. The size and shape of 
the markers are not to scale. 


pc. Assuming a flat universe, we compute the comoving dis¬ 
tance for each of the pointings, and the volume covered by 
the beam to that distance. The volumes corresponding to 
each extragalactic pointing are then summed to get the to¬ 
tal volume. With the above assumptions, the maximum red- 
shift probed in this survey was 0.17, and the average redshift 
among all our pointings is 0.13. We stress again the large as¬ 
sociated uncertainties. 

The average Galactic DM of all pointings that probed 
extragalactic space is 61.6 cm -3 pc, which, given our consid¬ 
erations about the host galaxy DM, justifies our DM spac¬ 
ing of 0.1 cm -3 pc (Figure pi. Given estimates of scatter¬ 


broadening related to these DM values in pulsars (Bliat 
|et al.|2004|), and com paratively small intergalactic scattering 
{L orimer et al.|20l3 1, we do not a priori expect scattering to 
render FRBs undetectable at the frequency and in the DM 
range of our survey. 

Assuming FRBs are standard candles with broadband 
emission, the volume sampled by the survey may depend on 
their spectral index rather than the maximum DM. This can 
be understood as follows. Sources at a given maximum red- 
shift should be bright enough at the frequency of observation 
to be detectable. Following Lorimer et al.]|2013]), we may be 
overestimating the volume in our survey if the spectral index 
of the sources is positive, as shown in Figure [5] where eq. 9 
from Lorimer et al. (2013) is plotted for various values of 
spectral index. Figure |5] shows the relationship of the peak 
flux density of a 5 ms FRB as a function of redshift and spec¬ 
tral index. Lines of different spectral index - 0.2, 0.1, —0.1, 
—0.2 - are plotted. A horizontal line is placed at the sensi¬ 
tivity limit of the ARTEMIS survey. For a given maximum 
redshift, the minimum spectral index of detectable FRBs is 
the one that leads to a curve that intersects the sensitivity 
line at that redshift. 

Figure [4] shows the sky coverage of the survey. Out of 
all the pointings in the survey, about 5% of pointings were 
limited to within the Galaxy due to the DM limit. We have 
excluded these pointings, leading to a total sky coverage in 



Figure 5. The relationship between the observed peak flux den¬ 
sity for 5 ms-duration pulses and the source redshift, based on 
treating FRBs as standard candles ( [Lorimer et al.|2013[ ) , for var¬ 
ious values of spectral index a. The black horizontal dashed line 
is the sensitivity limit of the ARTEMIS survey for 5 ms-duration 
bursts. The vertical lines are at the same redshifts as the vertical 
dotted lines in Figure [6] and help provide limits on the estimated 
spectral index of a = 0.1, as described in detail in the text. 
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Figure 6. Reproduction of the centre panel of Figure 2 from 
| Lorimer et al.| ( |2Q13| ). The white cross is the |Thornton et al.| ( |2013| > 
rate and the red arrow is the ARTEMIS upper limit. The white 
curve represents the FRB event rate at all redshifts constrained 
by the |Thornton et al7| ( |2013| ) rate, assuming that their survey 
was complete up to 2 = 0.75. The band around it represents the 
Poisson error on the |Thornton et al.| ( |2013) rate. The horizontal 
dotted line is at the level of the ARTEMIS upper limit, and inter¬ 
sects the extrapolated Thornton rate (and errors) at the vertical 
lines, which correspond to the redshift limits shown in Figure [5] 


our survey of 4193 sq. deg. The total duration of the survey 
as it pertains to sensitivity towards extragalactic bursts re¬ 
duces to 1446 hours. The single horizontal sensitivity line in 
Figure[5]is justified by the majority of our pointings being at 
high Galactic latitudes (>20°), where the sky temperature 
does not vary significantly. 
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3.1 Results 


Our 1446-hour survey resulted in no detection of FRBs. Fig¬ 
ure [7] shows an example of recorded events during the ob¬ 
serving session of 19 September 2014. These consist of sin¬ 
gle pulses from known pulsars and sources of RFI. During 
the survey, the real-time processing software detected single 
pulses from all expected bright pulsars, including B0950+08, 
B1919+21, B1133+16, B0525+21, B0531+21 (the Crab pul¬ 
sar), B2016+28, B2020+28, and B1237+25. This demon¬ 
strated the functionality of the system, including the abil¬ 
ity to process a large number of DMs in real-time. We 
also serendipitously detected signals from HAMSAT, a low- 
Earth-orbit satellite whose signals appear at low, non-zero 
DMs. Despite the sophisticated adaptive RFI excision algo¬ 
rithm, a small number of spurious events were recorded per 
epoch, which were all evaluated by eye. For a 5 ms resolution, 
the total number of such events accounts for approximately 
10” 4 of the total searched parameter space in DM and time, 
and is therefore negligible in all calculations. 

Based on the sky coverage and the time surveyed, and 
assuming an isotropic distribution of sources and complete¬ 


ness of the survey (for caveats, see Burke-Spolaor & Bannis- 
|ter||2014| |Keane fe Petroif||2015[ ), we derive an upper limit 
to the rate of FRBs at 145 MHz, of 29 sky -1 day -1 for 
5-ms-duration events above 62 Jy. 

The non-detection of FRBs at low frequencies and their 
low numbers at 1.4 GHz makes inferring the nature of these 
sources difficult. Given the estimated survey volume, we 
could have observed FRBs from 1.15xlO f) host galaxies, as¬ 
suming a galaxy density of 10~ 2 Mpc” 3 . The amount of time 
however that each of these galaxies was covered by the sur¬ 
vey beams is ~ 10” 3 yr, given the size and declination of 
the b eams. Assuming the FRB volumetric rate in Thornton 
et al. (2013) of 10” 3 Gal -1 yr” 1 , that leaves us with a single 


potential FRB event in our survey. We therefore have to con¬ 
sider the possibility that no FRB occurred during the survey. 
Given the uncertainty in the relationship between DM and 
redshift, and given the uncertainty in the volumetric rates 
of FRBs at the present day, it is useful to consider other 
possible reasons for the non-detection, such as the spectral 
index of FRB emission. 

As shown in Figure [5] where eq. 9 of |Lorimer et al.j 


(2013) is plotted for various values of the spectral index 


a, and as mentioned in the previous section, spectral index 
affects the volume surveyed. For values of a less than -0.2, 
the volume surveyed is dictated by the DM limit. For these 
values of a, following the standard candle model, we expect 
to see all bursts up to 2 « 0.13. 

Fi gure [6| r eproduces the centre panel of Figure 2 of 
Lorimer et al. (2013) up to 2 = 1, where the sky rate of 


bursts enclosed within a given volume is shown. The |Thorn-| 
ton et al. (20131 rate is marked by the white cross. Our 
upper limit on the rate, of 29 bursts sky -1 day” 1 , sits be¬ 
low the curve, and is marked by the red arrow. Assuming 
that FRBs occurred during our survey and given the caveats 
mentioned above, our null result could be explained by the 
most distant FRBs being too weak to observe, due to their 
spectrum. This leads to the realization that the maximum 
redshift is not set by DM considerations, but rather by the 
spectral index. To estimate the new maximum redshift, we 
observe that our sky rate is compatible with the |Thornton[ 


et al. 


(2013:) rate at a lower maximum redshift of 0.091 q'q 2 . 
TEis" is the redshift at which the curve of the scaled Thorn¬ 
ton rate of Figure[6]also becomes 29 sky” 1 day” 1 . FigureJH] 
allows us to associate a spectral index to the new maximum 
redshift probed, and errors on that spectral index based on 
the errors on the maximum redshift. Reading off Figure [5] 
this line of thinking allows us to place a lower limit on the 
spectral index of O.IIq'J. 

We must point out that this line of reasoning also re¬ 
duces the effective volume of the survey to 3.3 x 10 7 Mpc 3 , 
leaving approximately one third of the host galaxies we 
started off with that could host an FRB. If, however, the 
true reason for our non-detection lies in the spectral proper¬ 
ties rather than the volumetric rate, the implication is that 
the FRB spectral index is distinct from that of most pulsars. 
In the case of pulsars, a steep (negative) spectral index is 
thought to be indicative of a coherent emission process. The 
lack of such a steep spectrum in FRBs would serve as a clue 
to the nature of their emission process, which could be in¬ 
trinsically more narrow-band than pulsars, or characterized 
by a frequency dependent emission geometry. It is worth not¬ 
ing that two of the published FRBs at this time have been 
observed to have a positive spectral index within the band of 
observation: a most likely value of a between 7 and 11, con¬ 
sidering the frequency dependent telescope beam pattern, in 


Spitler et al. (2014), and 1.3±1.6 in Ravi, Shannon, & Jame- 


120151. Estimating spectral indices within an observing 
band is difficult due to effects of scintillation and frequency 
dependent gain variation within the telescope beam. This 
latter effect is more likely to render the observed spectrum 
to have a more negative slope than its intrinsic slope. There¬ 
fore, despite the small population of FRBs with a positive 
spectral index, a consistent picture may be emerging. 

In estimating the FRB spectral index from our survey 
we have assumed that FRBs are extragalactic standard can¬ 
dles with broadband radio emission. We have also implicitly 
assumed that the volumetric rate of FRBs, that reflects the 
rate per galaxy per year, is constant with redshift. We have 
not tested the standard candle hypothesis, while the vol¬ 
umetric rate may be a function of redshift, similar to the 
redshift dependent star formation rate. 


4 DISCUSSION AND CONCLUSIONS 

We have developed ARTEMIS, a GPU-based real-time fast 
transient detection system that works with individual LO- 
FAR stations. We have used the deployment at the Rawlings 
Array to demonstrate the functionality and performance of 
the system using known pulsars. We conducted a 1446 hour 
drift-scan survey for FRBs, searching the DM space up to 
320 cm” 3 pc. With a series of arguments relating to sensi¬ 
tivity and the intrinsic spectra of FRBs, and some simple 
assumptions about their nature, we estimate that we have 
surveyed a volume of 3.3 x 10' Mpc 3 . 

No FRB was detected in the survey at the Rawlings 
Array or indeed the Nangay station, and we derive an up¬ 
per limit of 29 FRBs sky” 1 day” 1 within the parameters of 
our search. The non-detection could be due to one or more 
of the following reasons: Not enough volume surveyed for 
enough time, interstellar and intergalactic scattering, and 
the intrinsic spectra of FRBs. 
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Figure 7. A sample output plot of single pulse detections in the ARTEMIS survey, showing a few pulsars and some RFI. Each cluster 
of markers correspond to single pulses from a pulsar, while the sets of points, elongated in DM, are RFI. 


Assuming the latter, we have attempted to constrain the 
spectral index of FRBs, to a > O.llg'j, assuming a standard 
candle model. The physical significance of this result is mod¬ 
erated by the fact that we cannot rule out other causes of 
non-detection at the frequency of our survey, including the 
possibility that FRBs are local events to the telescopes that 
have detected them. Assuming that the non-detection is in¬ 
deed due to the intrinsic shape of FRB spectra, our result 
is constrained by the sensitivity of the Rawlings Array. 

This work has implications for future low frequency 
FRB surveys. Firstly, at frequencies around 145 MHz, longer 
duration surveys with increased telescope sensitivity could 
allow further constraints of the spectral index to even higher 
values, or detect FRBs. It is however quite difficult to sam¬ 
ple substantially more volume of Universe, without running 
into greater potential scattering effects. Secondly, moving 
up in frequency by a factor of 2 and increasing the sensi¬ 
tivity by two orders of magnitude, which will ultimately be 
interesting for the low frequency component of the Square 
Kilometre Array, could allow surveys to probe redshifts up¬ 
wards of z= 1. Assuming scatter broadening is mainly due 
to propagation in the host galaxy and our galaxy, there is 
a reasonable expectation that it will not affect these predic¬ 
tions. This suggests that future high-sensitivity surveys at 
300 MHz, should result in a reasonable yield of new FRBs 
and should resolve the spectral index question. 

For ARTEMIS, we are developing a commensal real¬ 
time system in collaboration with the Berkeley SETI Re¬ 
search Center, for the Arecibo radio telescope. Named AL- 
FABURST, this system will use the seven-beam L-band feed 
of the Arecibo telescope to perform surveys for FR.Bs, po¬ 
tentially increasing the sample size. In this system, we aim 
to incorporate multi-beam coincidence rejection that will re¬ 
duce the number of false alarms, and automated triggering 
that will alert low-frequency telescopes in the event that an 
FRB is detected. Similar systems are being prepared for the 


South African SKA pathfinders (KAT-7 and MeerKAT) and 
for tied-array beams of the LOFAR International Telescope 
that offer an order of magnitude improvement in sensitivity. 
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